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Abstract. This paper presents the natural extension of Buckley-Feuring method proposed in [2] for solving 
fuzzy partial differential equations (FPDE) in a non-polynomial relation, such as the operator ip(D xl , D X2 ), 
which maps to the quotient between both partials. The new assumptions and conditions proceedings from this 
consideration are given in this document. 

Introduction 

Many approaches for obtaining non-numerical solutions of fuzzy differential equations (FDE) have been devel- 
oped from the introduction of fuzzy set concept by Zadeh [7\ . These ones give a diversity of definitions for FDE 
solution based on different notions of fuzzy derivative, such as Seikkala derivative, Buckley-Feuring derivative, 
Puri-Ralescu derivative, Kandel-Friedman-Ming derivative, Goetschel-Voxman derivative, or Dubois-Prade de- 
rivative . Some relations between these derivatives are presented by Buckley and Feuring in [3] . However, 
only a few of these fuzzy derivatives are valid in some contexts as FDE solution. For example, the Goetschel- 
Voxman derivative, or the Dubois-Prade derivative provide solutions that cannot be a fuzzy number, whereas 
Puri-Ralescu derivative, and Kandel-Friedman-Ming derivative, always exist and provide a fuzzy number as 
solution of the FDE, but making use of abstract subtractions of fuzzy concepts in their definitions, making 
difficult the interpretation of this solutions in some real applications. 

This paper uses the Buckley-Feuring derivative for solving FPDE. This derivative does not always exist, but if 
it does, provides a fuzzy number solution easily understandable in the context in which a specific FPDE has 
been developed. 

The authors which proposed this concept of derivative, developed a methodology for solving constant coefficients 
polynomial FPDE in [2] . This paper present the extension of this methodology to a non-polynomial expression 
in partial fuzzy derivatives. 



In the following lines, the components of a FPDE are enumerated: 

• Xi, i = 1,2, X\ 6 S\ C I\ = (0, Mi], X2 € S2 C h = (0, M 2 ]. Other domain limits can be 
established in this subsets, such as x\ > xi- 

• (3 = {Pi, /3a, /3fe), a triangular fuzzy number vector. 

• ^(Pj) is the membership function of (3j G j3j. 

• M/3j (a) = {(3j I f-i(Pj) > a, a G (0, 1)} set called a-cut. 

These sets are closed and bounded, so that is possible to define, for a fuzzy number j3j : $j[a] — 
[bi(a),b 2 {a)], where: 

— 61(a) is the lower value 0j in which fi(j3j) > a, f3j G 0j. 

— 62(a) is the higher value (3j in which yu(/9j) > a, (3j G [3j. 
and f3[a] = JJj 0j[a] 

• V(xi, X2, /3) is a positive and continuous function in (xi,X2) G Si x £2 with partials D Xl ,D X2 . This 
function must be also strictly increasing or strictly decreasing in x 2 G Si, that is V(k, xi, (3) is strictly 
increasing or strictly decreasing for all constant k G K. 

The fuzzy character of V{xi, X2, $) shown by the tilde placed over V is fixed by /3, and support the use 
of Buckley-Feuring derivative for solving FPDE. 
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• ip(D Xl , D X2 ) is an expression with constant coefficients in (D Xl , D X2 ) applied to V(xi, x 2 , (3). 

• F(xi,x 2 , (3) continuous function in (sc 1,0:2) G Si x S 2 . 

The specific FPDE treated in this paper has the following form according with this notation: 

<p(D Xl ,D X2 )V(x u x 2 ,P) = d V dXl = F(xi,X2,P) 

0V/OX2 

1. The Buckley-Feuring (B-F) solution 
The Buckley-Feuring (B-F) solution uses a solution of the crisp partial differential equation: 

V(xi,x 2 ) = G(xi,x 2 ,P), 

with G continuous y(xi,x 2 ) £ S± x S 2 . 
The next step is the fuzzification of G: 

Y(xi,x 2 ) = G(xx,x 2 ,f3), 

with G continuous V^i,^) £ Si x S 2 and strictly monotone for x 2 G S 2 . Note that Yi is only the fuzzy 
representation of G, but not necessary the solution of the fuzzy partial differential equation. If it finally 
happens and Y(x\,x 2 ) is a B-F solution, then V(x±, x 2 , f3) = Y(xi,x 2 ). 
With this notation, it is possible to see that: 

Y(x 1 ,x 2 )[a] = [yi(xi,x 2 ,a),y 2 (x 1 ,x 2 ,a)], and 

F(xi,x 2 ,(3)[a] = [fi(xi,x 2 ,a), f 2 (xi,x 2 ,a)],Va. 

and, by definition: 

yi(x 1 ,x 2 ,a) = mm{G(xi,x 2 ,f3), (3 e (3[a]}, 
y 2 (xi, x 2 , a) — max{G(xi, x 2 , (3), (3 £ 0[a]} ,and 

fi(xi,x 2 ,a) = mm{F(xi,x 2 ,f3), (3 E /§[«]}, 
f 2 (xi,x 2 ,a) = max{F(xi,a;2,/3), (3 G f3[a]}, 
Mxt,x 2 ,a. 

If it is possible to apply the f(D Xl , D X2 ) operator to yi,i = 1,2, getting continuous expressions V(xi,a;2) G 
Si x S 2 ,\/a, then it will be feasible to define the following expression in this domain r(xi, X 2 , a): 

T(xi,x 2 ,a) = [Ti(xi,x 2 ,a) 1 T 2 (xi 1 x 2 ,a)} 

with: 

ri(xi,x 2 ,a) = ip{D Xll D X2 )y 1 (xi,x 2 ,a) 
T 2 (x!,x 2 ,a) = if(D Xl ,D X2 )y 1 (x 1 ,x 2 ,a) 



For a B-F solution Y, it is necessary to be a fuzzy number for this one. If, for each pair (xi,x 2 ) G Si x S 2 , 
T(xi,x 2 , a) defines an a-cut of a fuzzy number, then, Buckley and Feuring [5] call to Yi(xi, x 2 ) differentiable, 
and we can write: 

ip(D Xl ,D X2 )Y(xi,x 2 )[a] = T(a;i,a;2,a), 

y{xi,x 2 ) eSiX 5*2, Va. 



So that, it is necessary to test if r(a;i, x 2 , a) really define an a-cut of a fuzzy number and verify the differentia- 
bility of Y(xi, x 2 ). For a triangular fuzzy number, the conditions are [1] : 

(1) tp{D Xl , D X2 )yi{x\, x 2l a) is an increasing function of a, for each (xi,X2) G Si x S 2 . 

(2) (p(D Xl , D X2 )y 2 (xi, x 2 , a) is a decreasing function of a, for each {xi,x 2 ) € Si x S 2 . 



THE EXTENSION OF BUCKLEY-FEURING SOLUTIONS FOR NON-POLYNOMIAL FUZZY PARTIAL DIFFERENTIAL EQUATIONS 



(3) ip(D Xl ,D X2 )y 1 (x 1 ,x 2 , 1) < (f(D Xl ,D X2 )y 2 (x 1 ,x 2 , 1) for each {xi,x 2 ) G Si x 5*2. 
Once delimited the differentiability concept of F(xi,:E2), it is possible to define the Buckley- Feuring solution. 
Y i (xi 1 x 2 ) is a Buckley- Feuring solution if the following conditions hold: 

(1) Y~i(xi,x 2 ) is differentiable. 

(2) ip(D Xl ,D X2 )Y(xi,x 2 ) = F(x u x 2 ,0). 

Obviously, if differentiability conditions hold by the candidate to B-F solution Y(xi, x 2 ), this one will be a fuzzy 
number. To complete the conditions for B-F solutions only is necesary to test that: 

<p(D Xl ,D X2 )Y(x 1 ,x 2 ) = F(x-i,x 2 ,0), 

or the equivalent condition: 

(1) ip(D Xl ,D X2 )y 1 (x 1 ,x 2 ,a) = F 1 (x 1 ,x 2 ,a). 

(2) ip(D Xl , D X2 )y 2 (xi, x 2 , a) = F 2 (xi,x 2 ,a). 
V(xi,x 2 ) 6 Si x S 2 , Va. 

In this case we can identify Y(xi,x 2 ) with V(x\, x 2 , /3). 

2. BOUNDARY CONDITIONS 

The FPDE can be subject to certain boundary conditions in a big variety of forms depending on a constant 
vector c = (ci, ...c„) e C\ x ... x C n . The inclusion of this ones has not great consequences in the methodology 
exposed. 

The crisp solution G acquires under boundary conditions the form G(x\, x 2 , (3, c). The fuzzification of G can 
take c in a triangular fuzzy vector c = (ci, ...c n ) E C\ x ... x C„ with c[a] — f] i A [a], i = l,...n and: 

Y(xi,x 2 ) = G(x 1 ,x 2 ,f3,c) 

In this environment with boundary conditions, it is necessary to add a new condition for a B-F solution: 
Y(xi,x 2 ) must satisficr these conditions. In this form Y i (xi 1 x 2 ) is a Buckley-Feuring solution if the following 
conditionshold: 

(1) Yi(xi,x 2 ) is differentiable. 

(2) <p(D Xl ,D X2 )Y( Xl ,x 2 ) = F(x u x2,0). 

(3) Y(xi,x 2 ) satisfies the boundary conditions. 

3. Example 

The following FPDE is proposed: 

%H7T =^i 1 x 2 , pG (0,1), m > l,x 2 >0 
ov jax 2 

In this example: 

F(xi,x 2 ,/3) = (3x^X2, e (0, l),xi >l,x 2 >0 
and, by definition, the operator 

tp(D Xl ,D X2 )V(x u x 2 ) — ► , 

A possible solution to this FPDE in a crisp environment is, without special boundary conditions: 

G(xi,x 2 ,/3) = x1x 2 +j,(3 e (0, l),xx > l,x 2 > 

with 0= (0,7) 
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Applying the fuzzification in j3 and 7, these ones acquire a triangular fuzzy number form and (3 = {(3, 7), (3 G 
(0, 1). While G holds all the conditions required, the Buckley-Feuring solution candidate is: 

Y(x 1 ,x 2 ) = G(x 1 ,x 2 ,(3,j) = xfx 2 +j,f3e (0, l),xi > l,x 2 > 

The fuzzy parameters have a membership function associated /i(/3) and ^(7) respectiveness. From the a-cuts, 
it is possible to define: 

(3[a] = [h (a), 6 2 (a)], 

7 N = 52(a)], and 

0[a] = (3{a] x j[a] 

And from this ones: 

Y(x-i,x 2 )[a] = [y 1 (x 1 ,x 2 ,a),y 2 (x 1 ,x 2 ,a)}, and 
F(xi,x 2 ,0)[a] = [f 1 (x 1 ,x 2 ,a),f 2 (x 1 ,x 2 ,a)],Va. 

Where: 

yi(xi,x 2 ,a) = mm{G(xi,X2,P), f3 G /3[a]} = G(xi, x 2 , gi(a), 61(a)), 
V2{xi,x 2 , a) = max{G(x 1 ,x 2 ,fi), f3 G /3[a]} = G(a;i,a:2, 52(a), 62(a)) 

and, 

/i(xi,x 2 ,a) = mm{F(x 1 ,x 2 ,(3), (3 G /3[a]} = F(x 1 ,x 2 ,g 1 (a),b 1 (a)), 
f 2 (x\,x 2 ,a) = max.{F{x\,xi,P), (3 G /3[a]} = F(x i; x 2 , 32(a), 62(a)), 

Va xi > 1, x 2 > 



In the G function proposed, \/{x\,x 2 )x\ >l,x 2 > 0: 

2/i(xi,x 2 ,a) = x b 1 l(a) x 2 + 51(a), 61(a) G (0, 1) 
2/2(2:1, £2, a) = x5 2(a) a;2 + 52(a), 6 2 (a) G (0, 1) 

f\(x- l ,x 2 ,a) = 61 (a)x^ 1 x 2 , 61(a) G (0, 1) 
h{xi,x 2 ,a) = b 2 (a)x^ 1 x 2 , 6 2 (a) G (0, 1) 

Testing the differentiability of G(xi, x 2 ,f3), from T(x\, x 2 ,a) — [ip(D Xl , D X2 )y\ (x\ , a; 2 , a), <f>{D Xl , D X2 )y 2 (xi, x 2 , a)] 
verifying if r(xi, x 2 , a) defines an a-cut of a triangular fuzzy number for each pair (xi, x 2 ), x\ > l,x 2 > 0, the 
following conditions must be hold: 

(1) tp(D Xl , D X2 )yi(xi,x 2 ,a) = b\{a)x^ l x 2 , 61(a) G (0,1) is an increasing function of a, for each pair 
(xi,x 2 ),xi >l,x 2 > 0. 

It will happen if tp(D Xl , D X2 )y\(xi,x 2 ,a) has positive derivative on a: 

d a (bi(a)x^ 1 x 2 ) — b\{a)'x^ x 2 

While (3 is a triangular fuzzy number and 6i(a)is defined from its a — cuts, 61(a) satisfies the condition 
and is an increasing function with 61(a)' > 0, and x^ x x 2 > for x\ > l,x 2 > 0, d a (61 (a)^ 1 ^) > 
and the condition holds. 

(2) tp(D Xl ,D X2 )y 2 (xi,x 2 ,a) — b 2 (a)x^ 1 x 2 , 62(a) G (0,1) is a decreasing function of a, for each pair 
(xi,x 2 ),xi >l,x 2 > 0. 

Again, it will happen if ip(D Xl , D X2 )y 2 {x\, x 2 , a) has negative derivative on a: 

d a (b 2 (a)xi 1 x 2 ) = b 2 (a)'x^ 1 x 2 
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We can now use that (3 is a triangular fuzzy number and 62(a) is denned from its a — cuts, so that, 
in analogy, 62(a) satisfies the condition and is a decreasing function with 62(a)' < 0, and x^ 1 x 2 > for 
xi > l,x 2 > 0, d a (bi(a)xi 1 x 2 ) < and the condition holds. 

(3) ip(D Xl ,D X2 )y 1 (x 1 ,x 2 , 1) < f{D Xl ,D X2 )y 2 (x 1 ,x 2 , l)Va;i,a;2) G h x J 2 . 
In the example: 

6i(l)a; ] ; 1 a;2 < b 2 {2)x^x 2 



This condition is holds automatically if we realize that 61 (a) and 62 (a) are defined from a triangular 
fuzzy number and 6 2 (1) > 61 (1) 

At this point, it is possible to say that Y(x\, x 2 ) = G(x\, x 2 , (3, 7) = x^x 2 + 7, (3 € (0, 1), x > 1 is differentiable 
and good candidate for B-F solution. But it is necessary the last step, and it must satisfier that: 

ip(D Xl ,D X2 )Y(xi,x 2 ) = F(xi,x 2 ,0) /3,7 € (0, l),V(x 1 ,x 2 ),x 1 >l,x 2 > 0. 
For G{x 1 ,x 2 ,f3,l) = xfx 2 +7,/3 e (0,1), a; > l,x 2 > 0: 

.DijY'Ca:!, 2:2) - /^xf 1 ^, (J G (0,1), an > l,ar 2 >0 

and, 

F(xi,x 2 ,0) = 0x^x 2 , G (0, l),xi >l,x 2 > 0. 

Thus G(xi , a: 2 , /3, 7) = x^x 2 + 7, /3 G (0, 1), x > 1, x 2 > is a Buckley-Feuring solution for this non-polynomial 
fuzzy partial differential equation. 
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